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Inflammatory responses in eucaryotic ceUs are often associated with oscilla- 
tions in the nuclear-cytoplasmic translocation of the transcription factor NF- 
kB. In most laboratory realizations, the oscillations are triggered by a cytokine 
stimulus, like the tumor necrosis factor alpha, applied as a step change to a 
steady level. Here we use a mathematical model to show that an oscillatory 
external stimulus can synchronize the NF-kB oscillations into states where the 
ratios of the internal to external frequency are close to rational numbers. We 
predict a specific response diagram of the TNF-driven NF-kB system which ex- 
hibits bands of synchronization known as "Arnold tongues". Our model also 
suggests that when the amplitude of the external stimulus exceeds a certain 
threshold there is the possibility of coexistence of multiple different synchro- 
nized states and eventually chaotic dynamics of the nuclear NF-kB concentra- 
tion. This could be used as a way of externally controlling immune response, 
DNA repair and apoptotic pathways. 

The synchronization between two oscillating signals exhibits a surprisingly deep level of 
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complexity [HI. Already in 1876, the dutch physicist Huygens observed that two clocks hang- 
ing on the wall tend to move in parallel after some time, i.e., they become synchronized [|2l. 
Since then, such phenomena have been observed in a variety of systems ranging from fluids 
to quantum mechanical devices IH |71 |6l HI El. In recent years it has become increasingly 
clear that living organisms offer a bewildering fauna of oscillators, e.g. cell cycles [|9l, cir- 
cadian rhythms [[TOl . embryo segmentation clocks [[TTIl . calcium oscillations [[T2l . pace maker 
cells [[T3ll . protein responses [fT4l [TSl . hormone secretion [fT6l . and so on. A natural question 
therefore is: do oscillators in cells, organs and tissues tend to synchronize to each other or to 
external driving oscillations? 

Here, we investigate the possibility of controlling the frequency of ultradian cellular oscil- 
lators by synchronizing them to external oscillations. Two important ultradian oscillators in 
mammalian cells are triggered by external stresses. After DNA-damage, the tumor suppres- 
sor protein p53 has been observed to oscillate with a period of 4-5 hours [fTTl [T8l . Secondly, 
inflammatory stresses have been found to lead to oscillatory behavior in the transcription fac- 
tor NF-kB fll5i . Bulk and single cell measurements after treatment with tumor necrosis factor 
(TNF) show distinct and sharp oscillations with a time period of 2-3 hours [fT4llT5l . 

In these two cases, oscillations are caused by a feedback loop which incorporates the forma- 
tion of a complex between the transcription factor and an inhibitor (Mdm2 in the case of p53, 
and IkBa in the case of NF-kB). The complex formation induces an effective time delay through 
non-linear degradation which suffices to generate oscillations in the transcription factor and its 
inhibitor, out of phase with each other [fT9l . Those two feedback loops have been modeled both 
by applying explicit time delays [|20ll and by modelling the complex formation [19|. An elab- 
orate model with 26 variables (mRNAs, proteins, complexes, etc.) for the NF-kB system was 
first formulated in Ref. lfT4ll . Krishna et al. reduced this model to the core feedback loop by 
assuming that complexes were in equilibrium ||2T|| . Fig. lA shows a schematic representation 
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of the resulting model that consists of three coupled non-linear differential equations: 

= kmn[Ntot - Nn) - kiiJ— — — — (1) 

at Kj + I Kn + I\n 



dt 



ktK - Imlm (2) 



^ = kulm - a[IKK]a{Ntot - A^n) (3) 
at Kj + 1 

The triggering stimulus, e.g. TNF, acts by changing the level of active IkB kinase, [IKK]a, 
which phosphorylates IkB, resulting eventually in its degradation. This degradation rate is one 
of the parameters of the model and Ref. [EH used different constant values of this parameter 
to represent different steady levels of the TNF stimulus. Default parameter values are given in 
Table 1. With these values and choosing [IKK]a = 0.5 jj,M as in [211 , one obtains sustained 
oscillations with a frequency u = 1/0.9 hr^^. 

Here we wish to examine the effect of an oscillatory TNF stimulus on the system. The sim- 
plest possibility is to assume the degradation rate of IkB would oscillate identically to the TNF 
stimulus, and Ref. [221 showed that this could result in chaotic NF-kB oscillations. However, 
unlike steady levels of TNF, representing an oscillatory TNF signal by a similarly shaped os- 
cillatory behavior of the IkB degradation rate is unjustifiable - non-linear interaction between 
TNF and IKK could well cause complex changes in the shape of the external signal as it is 
transduced. Therefore, we extended the model of Krishna et al [[2T1 to include the circuit that 
transduces the TNF signal to the IKK concentration, as shown schematically in Fig. IB. Ashall 
et al [[23][ have modelled this circuit in detail, and we add the two relevant differential equations 
from their model to the Krishna et al model: 



d[IKKl ^ k,^TNF]{[IKKU - [IKK], - [IKK],) - k,[IKK], (4) 



dt 



^^^^ = h[IKK], - kJiKKl , (5) 

dt ^"kA2Q + [A2Q][TNF] 



This model assumes that there is a constant pool of IKK which is interconverted between 
different states - active, inactive and neutral. TNF increases the rate at which inactive IKK 
is made active, and it is only the active IKK which phosphorylates IkB and thereby affects 
the degradation rate of IkBa. Ref. Il23l constructed these equations so that the TNF signal 
can be represented by a dimensionless number between (off) and 1 (on). We have used the 
parameters values from [|23l (see Table 1) except for Ntot and [IKK]tot- Ref. [23] takes both 
~ 0.1/iM, whereas [|2T| uses Ntot = 1/iM. Here, we chose to keep Ntot = 1/uM and varied 
[IKK]tot around l/^M, with \TNF] fixed at 0.5, to find a value that gave sustained spiky 
oscillations with a frequency in the range 0.3-1 hr~^ . The model in [|23l also includes another 
slow feedback via the molecule A20 as seen in the equations. For simplicity we ignore this 
feedback by keeping its concentration, [A20], constant as this feedback loop mainly fine-tunes 
the shape of NF-kB response and "there is a range of constitutive A20 expression values that 
can functionally replace A20 negative feedback" Il24ll . As default values, we finally chose the 
combination of [IKK]tot = 2 /iM and [^420] = 0.0026 yuM. This results in sustained spiky 
oscillations of frequency uq ^ 1/1.8 hr^^ when \TNF] is fixed at 0.5. Below, we show that 
with these parameter values oscillatory TNF stimuli in the model defined by equations (1)- 
(5) can produce both very organized responses - a multitude of synchronized states - as well as 
chaotic behavior. We have checked that our results are not qualitatively changed by perturbation 
of parameters around these default values, as long as the parameters result in sustained spiky 
oscillations when TNF is kept fixed at 0.5 (data not shown). 

We examine both sinusoidal as well as square wave oscillation of TNF, and in each case vary 
the amplitude and frequency of the applied stimulus, while keeping the rest of the parameters of 
the model fixed. When TNF is sinusoidally varied with different frequencies around the average 
value 0.5, we observe that the NF-kB oscillation synchronizes in interesting ways to the applied 
TNF frequency. For example, when the applied frequency is close to the uq, then the NF- 
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kB oscillation synchronizes to have the same frequency as the external stimulus, as well as a 
constant phase difference between the maxima of TNF and NF-kB (see Fig 2a). Similarly, when 
the applied frequency is in a band close to 2uq, the NF-kB oscillation synchronizes to a period 
exactly twice the applied frequency (see Fig. 2b). A fundamental result of our investigation 
is that the NF-kB oscillations will stay completely synchronized even if the frequency of TNF 
oscillations is slightly diminished or slightly increased! . That is, the external TNF signal is 
able to 'pull' the frequency of the NF-kB oscillation towards a rational ratio with respect to 
the applied frequency. This is known as phase (or mode) locking (31. As the amplitude of the 
applied oscillation increases, these bands of synchronization expand and the resulting shapes 
are called Arnold tongues jSl. Fig. 3a shows the Arnold tongues for the case where the TNF 
signal is sinusoidal, \TNF]{t) = 0.5 + Asin{2'Kvt), with varying frequency, v, and amplitude, 
A. Fig. 3b shows a similar diagram for the case when TNF is a square wave: [TNF] (t) = 
0.5 + A[sign{sin{2TTi't) + l]/2. This protocol may be easier to realize experimentally than a 
sinusoidal one, and in fact shows broader Arnold tongues. 

In principle, there is an Arnold tongue wherein NF-kB shows p peaks for every q peaks of 
TNF, for every rational number p/q (where p and q are natural numbers). The width of each 
tongue starts off infinitely small when A = 0, and expands smoothly as A increases. Evidently, 
this cannot happen without tongues overlapping, and indeed such overlaps occur as soon as 
y4 > 0. In general, in overlapping regions one expects to observe multistability, i.e. multiple 
synchronized states with different p/q values will coexist, with different states being realized 
when different initial conditions are used. However, for small A, the states corresponding to 
small p and q numbers generally dominate the observed behavior. As A is increased, these 
dominant states, such as 1/1 and 2/1, also start overlapping and then one can actually observe 
multistability (see Figs. 2c,d). As A is increased further, and there are more and more overlaps, 
one can also encounter chaotic behavior as shown in Figs. 2e,f. The same behavior occurs for 
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square wave oscillations of TNF (data not shown). 

This complex behavior of the existence, growth and overlapping of Arnold tongues is ob- 
served in several very simple sets of nonlinear differential equations, such as circle maps and 
other return maps (we refer the reader to E HI [51 for details). It has also been observed in a 
number of physical systems ranging from turbulent fluids, where synchronized states with ratio- 
nal numbers up to 83/79 haven been measured IH, quantum mechanical devices like Josephson 
junctions and semi-conductors [|26l [8l[6ll25l. crystals [7], and sliding charge-density waves (H. 
Synchronization is known to occur in living systems, such as fireflies, and circadian clocks en- 
train to the day-night cycle [|27l . However, to our knowledge such Arnold tongues have not 
been observed in vivo at a subcellular level. 

Our work suggests that this kind of intricate synchronization could be observed in the NF- 
kB system, and also the p53-Mdm2 system as it has a very similar core feedback loop. More 
specifically, we predict that: (a) oscillations in NF-kB can be synchronized to TNF oscillations, 
(b) the bigger the amplitude, the stronger the synchronization (when amplitudes are relatively 
small), (c) the oscillations can in principle be synchronized to all rational ratios with respect 
to the applied frequency, but states with smaller p and q values will dominate, (d) if oscilla- 
tions can be sustained for around a day, the states 1/2, 1/1 and 2/1 should be observable in 
practice, (e) when the amplitude of TNF oscillations is increased further, chaotic behavior will 
appear. Similar predictions could be made for the p53-Mdm2 system, controlled by external 
stimuli such as irradiation or exposure to DNA-damaging chemicals. If this basic synchroniza- 
tion technique works, we hope it could be developed into a tool to control inflammatory and 
apoptotic pathways and ultimately to regulate DNA repair, immune response and eventually 
cell fate. 
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Table 1 : Default values of parameters in the model. The first 9 are from Ref. L2JU and the next 4 
from Ref. [[23||. [IKK]tot and [A20] were chosen in order to obtain sustained spiky oscillations 
with frequency in the range 0.3-1 hr^^ when [TNF] is kept fixed at 0.5 (the actual frequency 
obtained with these values is z/q = 1/1.8 hr~^.) 
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Figure 1: (A) Schematic diagram of processes that form the core feedback loop controlling 
NF-kB oscillations in the model of [21 J. Nuclear NF-kB activates transcription of IkB which 
sequesters NF-kB in the cytoplasm. IkB kinase (IKK), when activated by external signals like 
the tumor necrosis factor alpha (TNF), causes eventual targeted degradation of IkB when it is 
bound to NFkB. The released NF-kB is then transported into the nucleus, closing the feedback 
loop. By assuming that all complexes are in equilibrium, these processes can be represented by 
three differential equations governing the dynahfiics of nuclear NFkB, IkB mRNA and cytoplas- 
mic IkB (see text). The external TNF signal's effect is simply modelled by appropriate choice 
of the degradation rate of IkB. This is sufficient when examining only steady levels of TNF, as 
is the case in fIT]. (B) However, when examining the effect of temporally varying TNF signals, 
this model has to be extended as shown to include the interaction between TNF and various 



Sinusoidal TNF stimulus 





Figure 2: (A) Simulation with applied TNF amplitude A = 0.09 and frequency v = 1/2.08 hr^^ 
(which is close to the reference frequency, = 1/1.8 hr^^). Resultant NF-kB oscillations 
have the same frequency as the applied TNF signal, giving the 1/1 "tongue" in Fig. |3j (B) 
Simulation with A = 0.09, u = 1/1.0 hr^^ (which is close to twice the reference frequency). 
Resultant NF-kB oscillations have twice the frequency of the applied TNF signal, giving the 
1/2 "tongue" in Fig. [3] (C) Simulation with A = 0.09, u = 1/0.83 hr^^ with initial conditions 
Nn = 1 /uM, Im = 0.5/iM and all others zdrfc fiM. The NF-kB frequency synchronizes to 
be 1/3 of the applied frequency. (D) Simulation with A = 0.09, u = 1/0.83 hr^^ identical to 
(C) except with different initial conditions where Nn = 0.1 fiM, Im = 0.5/iM and all others 
zero fiM. With these initial conditions, the NF-kB frequency synchronizes to 1/2 of the applied 
frequency. Thus, for this value of A and v, we observe multiple stable synchronized states. (E) 



(A,B) 




(unperturbed NFkB freq)/(TNF frequency) (unperturbed NFkB freq)/(TNF frequency) 



Figure 3: (A) Colors show the ratio of observed NF-kB frequency to the applied TNF frequency, 
as a function of TNF frequency and amplitude. The applied TNF signal is sinusoidal: [TNF] = 
0.5 + Asin(2Tiiyt). The Arnold tongues corresponding to the states 1/3,1/2,1/1,3/2,2/1,5/2,3/1 
are shown from left to right. For a tongue p/g, the boundary shown is the convex hull of points in 
the parameter space where the ratio of applied to observed frequency was within 0.5% oip/q. 
(B) Similar Arnold tongue diagram for the case where a square wave TNF signal is applied: 
[TNF] = 0.5 + A[sign{sin{2niyt)) + l]/2. 
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